Entanglement entropy of integer Quantum Hall states in polygonal domains 
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The entanglement entropy of the integer Quantum Hall states satisfies the area law for smooth 
domains with a vanishing topological term. In this paper we consider polygonal domains for which 
the area law acquires a constant term that only depends on the angles of the vertices and we give 
a general expression for it. We study also the dependence of the entanglement spectrum on the 
geometry and give it a simple physical interpretation. 
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I. INTRODUCTION 

The area law satisfied by the entanglement entropy of 
the low energy states of quantum many body systems in 
Condensed Matter and Field Theory, has become one of 
the most fundamental tools to study the physical prop- 
erties of these complex systems [i|-Q . To define this en- 
tropy one considers a low energy state f/;, usually the 
ground state, and computes the reduced density matrix 
PA by tracing out the degrees of freedom outside a do- 
main A. The entanglement entropy S'a('0)i associated 
to the state ip and the domain A, is defined as the von 
Neumann entropy of the reduced density matrix pA, i.e. 
Sa{^) = ^^T^PA^ogpA- The area law states that Sa{'4') 
is proportional to the size of the boundary of A. In 3D 
this size is the area separating A from its environment, 
which gives the law its name. In lower dimensions one 
should rather used the terms perimeter law in 2D and 
zeroth law in ID, but those names are not customary. 

Three issues are important regarding the area law: vi- 
olations, fluctuations and subleading corrections. They 
all provide a great amount of information about the sys- 
tem. In conformal invariant ID models, the area law 
shows a log violation proportional to the central charge 
of the corresponding CFT and the topology open/close 
of the system [4H7'| • Fluctuations around the log law for 
the Renyi entropy in Luttinger liquids allows a deter- 
mination of the Luttinger parameter [8| , and subleading 
terms contain information about the scaling dimensions 
of the operators [0] . 

In 2D much less is known about corrections to the area 
law. In systems with topological order, the correction 
is a constant term 7 denoted topological entanglement 
entropy [ifl, [U , and whose value is given by the loga- 
rithm of the total quantum dimension of the anyonic ex- 
citations. The area law of the Fractional Quantum Hall 
(FQH) states has been the target of several recent studies 
|12h16| |. in order to confirm its validity and to compute 
the value of 7 predicted in [l^, [ll| ■ Reference [lg| uses 
Chern-Simons theory, finding the predicted value of 7, 
however the linear behaviour of Sa, is not captured, due 
to the purely topological nature of this theory. There are 
numerical studies using the Laughlin wave function [12- 
Il4| and exact diagonalization ,1\ ; ISj , for filling fractions 
1/ = 1/3, 1/5 and the v = 5/2 Pfaffian state. The ap- 
proaches of [13 - [l5| use the orbital basis for the Landau 



levels. The close relationship of this basis to the spa- 
tial partioning of the blocks leads to an area law of the 
form Sa = c^/Ta — 7 + 0{1/Ia), where Ia is the num- 
ber of Landau orbitals in the block A. The numerical 



values of 7 computed in the spherical geometry [12| and 
the torus geometry [14l.ll5| agree, within some precision, 
with their theoretical values, despite of the fact that the 
systems analyzed are not very large. 

For non topological models there are also some results. 
In those with a Fermi surface, the area law exhibits a log 
violation |17H23 l reminiscent of the ID conformal sys- 
tems, which suggest that higher dimensional bosoniza- 
tion methods may give an appropiate description [23|. 
Corrections to the area law in non-smooth domains, in 
2-1-1 dimensions, have been obtained in relativistics free 
bosons and fermions in ref. [2J] and in an interacting CFT 
using the Ads/CFT conjecture in [2q]- The 2D quantum 
critical models of Fradkin and Moore [2a] with critical 
exponent z — 2 have also universal subleading contribu- 
tions (see also i27i|). 

The aim of this paper is to calculate the entanglement 
entropy of the IQHE ground state in arbitrary polygonal 
domains. We obtain that the entanglement entropy is 
given by the area (or perimeter) law plus a contribution 
due to the vertices, 7, that only depends on the angle of 
the vertices and the density of the fiuid of electrons. 

In the following we shall concentrate in the IQHE state 
with filling fraction ^ = 1 defined on a cylindrical geom- 
etry, but the calculation of 7 for higher filling fractions 
and different geometries are easily generalizable. 

Let us consider the Landau model for a particle in a 
cylinder of size L^xLy. The one particle wave function in 
the lowest Landau level (LLL), in the gauge A = B{0, x), 
is (in units of the magnetic length i equal to one): 



4>ky{X,V) 



1 



711/44/2 



pikyV „~{x-ky) jl 



(1) 



On the cylinder, the identification of the wave function 
along the y direction implies: 
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^y- L. 



"0,1^ , no 

h 1 < n < — . 

2 2 



(2) 



The number of LLLs, no, is obtained imposing that the 
particle lives in the strip \x\ < Lx/2, which yields uq — 
"" - . This value also gives the total number of quantum 
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Zl 



FIG. 1; Domain X> = Va U Vb including two vertices with 
angle a. Varying the parameter l^ ~ — °*g we can obtain 

the correction to the area law, 7, for a £ [0,7r]. 



fluxes through the box. The electron operator can be 
written as 



i^{x, y) = ^ ipky [x, y)cky + higher LLs, 



(3) 



where H{x) = — xloga; — (1 — a::)log(l — x). We have 
introduced a parameter l^ that allow us to vary the angle 
a (see fig[T|). If l^ = 0, the domain 'D{lx = 0) = I?o 
becomes a half cylinder. The eigenvalues of the correlator 
([7]) for this case were given in [29{ : 



A. 



•D( 



4'n{x,y)(l)n{x,y)dxdy = - 



I- Erf 



/2mT\ 



-00, .., 00 



(10) 



In this formula we assumed that Lx is effectively infinite. 
Computing now the entanglement entropy © using (jlOp 
one obtains St>o — 0.203291 * Ly. Observe that Ly co- 
incides with the perimeter of Vq, but this is a general 
result that holds for any smooth domain 'Dsmootin i-e.: 

5p_.. =c*P, (11) 

where P is the perimeter of T) smooth and the constant 



c = 0.203291 



fl2) 



is independent of the geometry of the system but varies 
with the number of fully occupied Landau levels [23] . 



where Ck is the fermionic destruction operator of the 
LLL labeled by ky. The extra term in ([3|) involves the 
remaining Landau levels, which are empty for filling frac- 
tion V = \. The ground state for z/ = 1 is given by: 



l'i'o)=n,„4 JO), 



(4) 



where |0) is the Fock vacuum. The two point fermion 
correlator in this state is, 

Cry^{<^o\i^\x,y)i^{x',y')\'^o). (5) 

Using ([3]) and (j4]) one finds: 

Gr,r'-5]</':,(x,2/)(/.fc,(x',2/')- (6) 

ky 

We want to compute the entanglement entropy, Sv, of 
the state $01 in a polygonal domain T) embedded in a 
cylinder of radius Ly such as that shown in figHJ This 
entropy is given by the formula Sx> = — tr px> log p-p where 
PT) — trx)^|<&o)(^o|i and Vc is the complement of I? in the 
cylinder. The computation of Sj} is done in two steps 
@i[2M|- First one restricts the correlation matrix Cr.r', to 
the domain T), i.e. 



Cr,r'=Cr,r', r, r' S 2?. 

Next, one diagonalizes Cr,r', i-G- 

d^r' Cr,r' 5(r') = A5(r). 
The entropy ^x; is obtained by means of. 



L 



(7) 



(8) 



(9) 



II. CORRECTION TO THE AREA LAW 

The aim of this section is to show that the area law for 
a polygonal domain T) is given by 



Sv = S-D,^„„th + 7, 



(13) 



where S-D^^^^^t^ is the entanglement entropy of a smooth 
domain Vsmooth with the same perimeter as P, given by 
eq. (fTTj) . and 7 is a constant term due to the corners of 
the domain. 

To proof eq.([T3]) we start by diagonalizing the correla- 
tor ([7]) in a generic domain V. Using ([1]) and © equation 
(|5]) can be written as: 

^(^;(a;,y) / 0„(x',2/')5(a^',2/') = ^ C(a^>2/)^n 

n ''^ n 

= A5(x,y) , (14) 

where we have defined: 



A„ 



dx'dy'(j)n{x',y')g{x\y') 



(15) 



■D 



The vanishing eigenvalues A of equation (|S]), does not 
contribute to Sv, so one can focus on the non vanishing 
ones. Using equation (1141) . the corresponding eigenfunc- 
tions g{x,y), can be written as 



9ix,y) = j^<j)*^{x,y)A„ 



(16) 



Finally, replacing P^ into (ITSt one obtains the following 
equation for the eigenvalues: 



I](^^)n-^'"=^^» ' 



(17) 
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FIG. 2: Different domains in tlie torus with 4,6,8 and 10 rigiit 
angles. 
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FIG. 3: Plot of a{a) vs. a for a G [0,7r]; the dotted line 
corresponds to the values of a{a) obtained by diagonalization 
of the Ft) matrix (|17|) and the continuous line shows fit 
There is a perfect matching between the two curves. 



where F-p is a matrix with elements 



iFv)r. 



4>n{x,y)(t>*n{x,y). 



(18) 



V 



Observe that the continuous eigenvalue equation ([5]) has 
been converted into a discrete one pT)) . where the eigen- 
values of the matrix F-d are those of the correlator ([7]) • 
This fact will allow us to apply numerical methods to 
find the eigenvalues A. 

Indeed, let us first show the validity of equation P^ 
for some simple domains as those depicted in figlH The 
domains I?„, with n = 4, 6, 8, 10, are four different poly- 
gons with n right angles. For each of these domains we 
first compute the matrix Fx> fsee llSp . and we diagonalize 
it for different values of P. The corresponding entropies 
Svn can be fitted to the formulas: 



a) Svi = 

b) Svf, = 

C) S-Ds = 
d) Sv^a = 



-0.341315 -I- 0.203291 * Pt>^, 
-0.512909 + 0.203295 * Pv^, 
-0.683902 + 0.203291 * Pv,, 
: -0.85412 + 0.203289 * Pv,o- 



Notice that the constant multiplying the perime- 
ters Pt)„ agrees, to great accuracy, with the value c in 
dill) for smooth domains. From these examples one can 
extract the value of 7 for right angles: 

7 ( — j = a ( x ) * (number of 7r/2 vertices) , (19) 

with a (1) = -0.0855. Based on (|13I19|) we can pro- 
pose a general expression for the entanglement entropy 
of an arbitrary domain V with n vertices parameterized 
by angles at {i = 1, . . . ,n) 



Sv 
7({«J) 



E 



■7, 



with 



a{ai) *n.i , a^ G [0,27r], 



with rii the number of vertices with angle a^ in the 
boundary of V. Each of these angles varies between 
and 27r, but as we show below, the parameter 7 is invari- 
ant under the symmetry a o 27r — a, which follows from 
the equality of the entropy of a pure state in a domain 
and its complement. For this reason we can restrict a to 
the interval [0, tt]. 

The calculation of a(a) can be done considering the 
domain of figUl For different values of a S [0,7r] (or l^ 
in figHJ, we diagonalize the F-d matrix, calculate Sv and 
finally using (j20p . we obtain the value of a{a) by means 
of the equation 



a(a) ^\{Sv- 0.20329 * Pv) 



(21) 



with Pt) the perimeter of T). The factor i in ([21]) arises 
from the fact that the domain of figH] contains two ver- 
tices with angle a. The numerical determination of a(a) 
is given in figlH This curve can be fitted with the follow- 
ing expansion in the variable 21x1 Ly — |cot(Q:/2)| (see 
figffl): 

a{a) = 0.00723706 |cot(a/2)| - 0.155709 |cot(a/2)|^ + 
0.0876488 |cot(a/2)|^ - 0.0306014 |cot(a/2)|'' + 
0.00698332 |cot(a/2)|^ - 0.00105786 |cot(a/2)|^ -I- 
0.000105322 |cot(a/2)f - 6.616 * 10"^ |cot(a/2)|* + 



2.375 * lO"'^ |cot(a/2)|^ - 3.7112 * lO'^ |cot(a/2)| 



10 



(20) 



(22) 

Note that ([^^ satisfies the relation a{a) ~ a(27r — a) that 
together with ((20)) imply that S-d — St>^ {T>c is the com- 
plement of T)) as explained above. Hence, we can restrict 
ourselves to the interval a S [0, tt] where cot(a/2) > so 
that we can drop the absolute values in ((22)) . 

In order to assess the accuracy of the fit ((22p we can 
compare the entanglement entropies obtained by diag- 
onalization of the Fj) matrix ((TT)). for the complicated 
domains of figJH with the theoretical formula ((^0)) using 



a) 
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FIG. 5: Overlap between the eigenvalues \n{lx ~ 21, Ly = 30) 
obtained by diagonalization of the Ft> matrix (|17p and those 
obtained using the deformation (|25p with k,{Ix, Ly) = 0.40266. 



ous results have been obtained for a cylindrical geometry, 
however it can be shown that a{a) is the same function 
for the sphere and the plane and thus 7 is also indepen- 
dent of the geometry. 



C) 




— t 1 1 1 1 1 1 — 

■30 -21 -14 -7 7 14 21 30 



FIG. 4: The domains a) b) and c) with perimeters P=202.479, 
260.45 and 254.3033 respectively are used to check the validity 
of (I20I) and (I22I) . 



the extrapolation ([2^ . In the case of diagonalization of 
the Fx) matrix we obtain: 

Sa ^ 38.04348652 , Sb = 50.02330873 and 

Sc = 48.92750022 (23) 

and using ((20|) and ^22\i we arrive at 

Sa = 0.20329* P + 4*a(7r/2) + 6 *a(arctan(9/15)) 
= 38.047, 

Sb = 0.20329*P-f 2*a(7^/2) + 4*a(arctan(9/15))^- 
5 * a(2 * arctan(9/15)) = 50.0253 and 

Sc = 0.20329*P + 2*a(7r/2) + 2*a(arctan(9/15))-f 
2 * a(arctan(7/15) + arctan(9/15)) + 
5 * a(2 * arctan(7/15)) = 48.927. (24) 

The comparison of (P5)l and ([M|) shows clearly the va- 
lidity of (pHI and of the extrapolation ([22). The previ- 



III. ANALYTICAL CALCULATION OF a{a). 

In this section we shall give an analytical justification 
of equation ^U\\ . A crucial point in our analysis is that 
the eigenvalue Xn{lx,Ly) of the Fx> matrix ([T7)) . corre- 
sponds to a deformation of the eigenvalue ^TU\\ for the 
case of the smooth domain T>o [l^ = in figll]). This 
result can be shown by diagonalization of _Fp for many 
different values of Ix and Ly. This deformation is param- 
eterized by a function K(lx,Ly) and it is given by : 



An(ia;, Ly) 



I- Erf 



2mr 
L11 



(1 - K{lx,Ly)) 



(25) 



In figE] we show, for the particular case Ix = 21 and 
Ly = 30 (a = 71.075° in figlS]), the perfect matching 
between the eigenvalues obtained by diagonalization of 
the Fx) matrix with those given by (|25t . 

To calculate the function K,{lx,Ly) we can use equa- 
tions ([5)l. (pni) . (|25l) and the change of variables x = 
(2Tm/Ly) (1 — K{lx,Ly)) obtaining: 



/oo 
H[\n{lx,Ly)]dn 
-CX3 



LyC 



(1 - K.{lx,Ly)) 



(26) 



with Pxi = 21/ ( ^ ) +1"^ the perimeter of V and c given 

in (fT2|) . From equation (p6| . the deformation K{lx,Ly) is 
expressed by: 



l^{l-Xl Ly) 



1 



LyC 



{cPv + 2a(a)) ' 



(27) 



with a{a) given by the extrapolation ((22|) or by the the- 
oretical value (15^ to be obtained below. 

Equation (|25|) is rather interesting since it shows that 
the whole spectrum of the matrix Fx> is given in terms, 
basically, of the function nil^, Ly). Hence, we expect that 
overall quantities like Tr F^ — ^nXni^x,Ly), will also 
depend on the latter function. If for some p, we were 
able to compute this trace, then we will know K{lx,Ly) 
as well a{a) by means of eq. (P7|) . The simplest choice is 
p = 1, but in this case the trace of Fx> does not depend 
on k{Ix, Ly). This follows from the relations 



TrFr, = ^ A„(?„ Ly) = / C{x, y)dxdy ~ pqA, (28) 

where C{x,y) is the correlator d?]), po is the density of 
electrons and A >> 1 is the area of the domain V. 

The next choice is p = 2 for which Tr F^ does depend 
on K{lx,Ly). Indeed, using ([25]) and taking the limit 
?io >> 1 we obtains: 



nO 



TrF%^ J2 ^l{l.,Ly) = 



n——nO 



no 



2V27:3/^-l + Kih,Ly)) 



(29) 



Then from ([291) and ^, the entropy Sv and TrF^, are 
related by 



-2V2Tr^^^c (TrFl - 



no) 



(30) 



Therefore if we find an alternative formula for TrF^, we 
would obtain a theoretical formula for S-p and justify in 
this way equation (|20p . This is done in the appendix for 
the domain T> of fig[T]with the result: 



TrF^ 



no 



2V27r3/2 
4 + (tt - a){-l + 3Cos{a))Csc{a) 



(31) 



where Pp is the perimeter of the domain V. Finally, from 
([50)1 and dnj) we obtain: 

Sv = cPv + 7(0;) with 
c 




a 



FIG. 6: Plot of a{a) vs. a. Continuous curve corresponds to 
the extrapolation of a[a) eQ. (|22|) obtained by exact diagonal- 
ization of the Ft> matrix. Dashed curve corresponds to the 
theoretical approximation for a{a) in ea. (|32p . 



IV. ENTANGLEMENT SPECTRUM 

The reduced density matrix V of a domain contains 
of course much more information that the entropy St>- 
In the context of the QHE this information is related to 
the physical degrees of freedom of the edge excitations, 
as proposed in reference [30| . This is a rather surprising 
conjecture since, after all, the edge of V is an arbitrary 
curve within the whole system. The entanglement Hamil- 



tonian, H-d, defined as pT> = 



— f^ — HT} 



is then expected to 



be intimately related to the Hamiltonian describing the 
excitations of a real edge. In the QHE the edge excita- 
tions are described by a chiral CFT, which also describes 
the excitations of the bulk. For the IQHE this result 
follows from the entanglement spectrum that is given by 
the eigenvalues of the operator [3l| 



(1/2) - C; 



■D: 



(33) 



with C-D the two point correlator defined in ([T]). In the 
smooth domain Do (obtained considering l^ — in figH]), 
the eigenvalues of Cx>o are given by eq. ([TU]) , so the entan- 
glement spectrum will be: 



-Erf 



27rn 



-00, 



., CX) 



(34) 



7(a) 



2a{a) = 
csc(a)) 



/2Tr 



with qn the eigenvalues of Qvo- In the limit Ly » 1, 
(4 + (tt — a)(—l + 3cos(q;))x this spectrum becomes g„ = 2Y^n/Lj,, which is that of 

a free chiral boson moving on a circle of length Ly with 
(32) a velocity v — I/a/tt- If we now deform the domain as in 
fig. [21 the entanglement spectrum in the boundary of V 
can be obtained from eq. (|25p 



Fig. ([6]) shows the overlap between the extrapolation a{a) 
in (1221) obtained by diagonalization of the F-p matrix and 
the theoretical expression ([32)) . The small difference be- 
tween the curves is due to the approximation done in the 
calculation oiTrF^ (see eq. ([45l) in the Appendix). 

In conclusion, in this section we have proof the valid- 
ity of our proposal equation (|20p for the entanglement 
entropy in a non-smooth domain. 



-Erf 



-00, . 



27rn(l — k{Ix, Ly) 

l:, 



(35) 



In the limit Ly » 1, eq.([17|) implies P-p = Ly/{1 
K(J,x,Ly)), and the spectrum ([55)) becomes g^ 



^^pKTijPx) which is that of a free chiral boson moving 
on a circle of length P-p with a velocity v = I/i/tt- That 
circle has the same length as the boundary of Pp. The 
latter boundary has two singular points on it, i.e. two 
vertices with angles a and 27r — a, but they have sub- 
leading effects both on the entropy S-d and the entan- 
glement spectrum. This phenomena is analogue to what 
happens for real systems where the electrons surround 
the hills and valleys of the potential travelling all the 
way through the sample. 



Appendix 

In this Appendix we calculate the trace (p9)) analyti- 
cally. Let us first write it as 

TrFl^uT,^ = TrFl^ + TrF^^ + 2Tr (Pp.PpJ , 

(38) 

with T^A and Db given in figHJ From the definition of 
Fv (see eq.lfTT])). the first term in (|38|) reads 



In summary, we have computed in this paper the sub- 
leading term to the area law for IQHE states on polygonal 
domains, which is given by the expression 



= ^a(aj) *ni 



(36) 



where n^ is the number of vertices in the domain with an- 
gle ai and a{ai) is a function that only depends on the 
angles of the vertices and the density of the fiuid. For the 
IQHE with filling fraction v = l,we have found numeri- 
cal and analytical expressions of the function a[a), given 
by eqs. ((22)) and ([32]) . which agree rather well. The fact 
that the correction ([36]) is a constant has its origin in the 
gapped character of the IQHE. For gapless 2-1-1 systems 
one may expect a size dependent subleading term. In- 
deed, this has been confirmed for relativistics free bosons 
and fermions in ref.[2J| and in an interacting CFT us- 
ing the AdS/CFT conjecture in [2^. Also, the critical 
models with dynamical exponent z = 2 of Fradkin and 
Moore [2a], exhibit a logarithmic subleading correction 
of the form 










ra{x' ,y')4>n{x' ,y') X 


0„(x,2/)0™(a;,y) ^^J J e-|(--')'-^('^-^')' = 


/■^ f g-^ix-x'f~^a^{l^+x+x'f g-i(4+a^)(x-x')^ 


l^'^"^^"^ l^ 2^2 27r2 


ae~i(^~'=">\x' ~x)Erf 


a{x — x') 


2 


2v^ 


ae 2(=^ ^) 

1 ^^ 


4^2^3/2 


' 4\/2^3/2 


{l,+x + x')Erf 


a{lx + X + x') 
2^2 


) 


=A+B+C+D , 



(39) 



where a — j^. The integrals A and B in ([5^ can be 
easily found 



A = 
B = 



7r-|-2arctan[i - f] 



An^a 



V^lx 



(4 + a2)7r2 ^74^^7^3/2 



(40) 



To compute the integral C we first integrate the variables 
y, y' and x obtaining that: 



/3 = /(a) * LogiL) 



(37) 



with L the perimeter of the domain and a the angle of 
the vertex. The function /(a) has a similar behavior as 
a{a), in the sense that both curves satisfy /(a), a{a) < 
and g^' , °^^^' > 0. These properties are a consequence 
of the strong subadditivity relation satisfied by the en- 
tanglement entropy [25j . 

Finally, we have obtained the entanglement spectrum 
of non smooth domains which corresponds to a chiral 
free boson moving on the boundary, in agreement with 
the conjecture of reference [3(l|. 
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C 



dx 



2V27r3/2 



,-\(l:,-2x'f 



Erf 



a{lx — 2x') 



4V2 



27r3/2^2(4-|-a2) 
C1+C2. 



Erf 



4^2 



(41) 



The integral Ci only depends on the ratio a = Ly/l^. 
Therefore if we vary l^ , maintaining the value of a con- 
stant (by adjusting Ly), the integral C'l remains the 
same. Taking the limit Z^; — >■ 00 in Ci and making the 
change of variables u = Ix — 2a;', one obtains an integral 
between and 00 whose value is 



Ci = — - arctan 

27r2 



The integral C2 is easily done 

a2 a^l^ 



C2 



(4-fa2)7r2 2^2(4 -Ha2)7r3/2 



(42) 



(43) 



Collecting the previous expressions one gets 



to obtain and they read 



C = ——^ arctan 



(4 + a2)7r2 2^2(4 + a2)^3/2 



(44) 



The integral D in ^: 

r-if ae-^^''-''">\l^ + x + x')Erf 



D = 



a{lj:+x-\-x ) 
2\/2 



4V27r3/2 



(45) 



cannot be solved analytically. However, one can obtain 
a very good approximation replacing in ((45]) {Ix + x + 



x')Erf 
tion yields 



a(la:+x+x') 
2\/2 



D = - 



by (Ix +x + x'), which upon Integra- 

(46) 



air 



all 



2V27r3/2 47r ' 
Finally, from (gOl) , dMl) and dH]) we arrive at 



V^lx 



TT + 2 arctan [\ - f ] 



47r2a 



\/4T^7r3/2 27r2 



arctan 
alx 



(4 + a2)7r2 
a2 
^ (4 + a2)7r2 " 

a?? 



2^2(4 + a2)7r3/2 2N/27r3/2 47r 



(47) 



To complete the calculation one needs the quantities 
TrF^^ and 2Tr {Fv^F-d^) in ([38]). Their values are easy 



TrF^s = E (^^b)„„ (^i^J„„ = "0 - ^^^3/, ~ ^ 



L 



L 



TrF(p^T>B) - Zl (^T^A)nn (Fv b) nn - ^"^f^j^ g/ j^ ' 
n V 2^ 

(48) 

with riQ defined in ([2]). In (|48l) wc have used the fact that 
in the domain T^b, the matrix F is diagonal. 

Finally, from (gT]) and (gH]) we obtain that: 



TrF, 



F, 



VaUVb 



a 



no 



VauVb 2y27r3/2 ' ■"" ' (4 + a2)^2 ' (4 + ^2)^2 

1 a arctan [i — |] a arctan [|] 



a 
4a7r 47r 



2a7r2 



27r2 



(49) 



with P-DauVb the perimeter of the domain V in figH] 
Observe from fig[T] that the ratio a 



-r- is related with 

ix 

the angle, a, of the vertex by a = 2Tan{a/2). Therefore 
we can write eq. (|49p in terms of the angle a obtaining 
that: 



FrF-p^^j,^ 



Pi 



VaUVb 



no 



2y27r3/2 
4+ (tt — «)(— 1 + 3cos(a))csc(a) 



(50) 
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